Brownian molecular motors driven by rotation-translation coupling 
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We investigated three models of Brownian motors which convert rotational diffusion into directed 
translational motion by switching on and off a potential. In the first model a spatially asymmetric 
potential generates directed translational motion by rectifying rotational diffusion. It behaves much 
like a conventional flashing ratchet. The second model utilizes both rotational diffusion and drift to 
generate translational motion without spatial asymmetry in the potential. This second model can 
be driven by a combination of a Brownian motor mechanism (diffusion driven) or by powerstroke 
(drift driven) depending on the chosen parameters. In the third model, elements of both the Brow- 
nian motor and powerstroke mechanisms are combined by switching between three distinct states. 
Relevance of the model to biological motor proteins is discussed. 

PACS numbers: 87.16.Nn, 87.17.Aa, 05.40.Jc 



I. INTRODUCTION 



The mechanisms of motor protein motility have been 
a major research topic in biophysics and biomedical re- 
search for several decades 1]. Pioneering work by Hux- 
ley and his co-workers |3, y| established a mathematical 
model for muscle fiber contractions based on a so-called 
powerstroke P, |3 involving the conformational change 
induced by ATP hydrolysis in the cross-bridges formed 
between the myosin and actin filaments. This conforma- 
tional change in the neck region of myosin has been con- 
firmed by various experiments 0, lu 13. H S Il3 ■ Based on 
Huxley's model, it has been believed that this conforma- 
tional change induces a s wing ing motion of the myosin 
(lever arm hypothesis) [l|, lll| , rotational motion about 
the neck which is in turn converted into linear transla- 
tional motion. This model suggests that a motor protein 
utilizing a powerstroke would behave like a stepping mo- 
tor with a fixed step length per ATP consumption. In 
addition, the step length is expected to be proportional 
to the neck length of the proteins. 

Recent single molecule experiments testing a variety of 
protein mutants with different neck lengths support the 
lever arm mechanism of protein motility |12. Il3l UA . llSl 
lla. llTl Il8l | . Other recent experiments have emphasized 
both the length of the motor protein neck as well as the 
magnitude of the angle through which it swings during a 
powerstroke in an attempt to explain the unusually long 
step length of Myosin Id [l^- However, some similar 
experiments showing unexpected step sizes for a given 
neck length appear to be inconsistent with the lever arm 
model ,20!, :2i |2l |2l IM ISl IM 123, IM l^l. In these 
experiments, proteins must somehow move over several 
steps during a single chemical cycle in order to produce 
the observed step size. Some experiments showed that 
the speed of the motor protein is independent of the neck 
length in contradiction with other experiments |23| . Fur- 
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thermore, the step length of certain motor proteins seems 
too long compared with their neck length 13(1 l31l and 
some diffusive steps appeared to be necessary [llll33.l33J . 

The filaments along which motor proteins move are 
polymers made up of polar units with a "plus" and a "mi- 
nus" end. Most myosin-based motor proteins move pref- 
erentially towards the plus end of the actin filament track. 
Myosin VI moves in the opposite direction 27, 28, 31, 34j. 
In addition, other motor proteins such as Myosin IXb re- 
cently have been observed to move in the opposite direc- 
tion under certain conditions 35, 36]. These anomalous 
motor proteins are structurally very similar to the their 
plus-ended counterparts. It has been observed that when 
the lever arm is engineered in a certain way, the mo- 
tor moved in the opposite direction |33|- In the case of 
Myosin VI, it has been speculated that the naturally oc- 
curring protein moves preferentially towards the minus 
end of the filament through a similar mechanism |34| . 
However, the actual mechanism determining the direc- 
tionality still remains unknown. 

The majority of motor proteins are comprised of two 
motor domains that work in conjunction with each other 
to produce directed motion. The two heads are partic- 
ularly crucial to the models put forth to describe pro- 
cessive motors such as Myosin V, which move through 
many ATP hy drolysis cycles before detaching from its 
filament [Ellg. Most of the prevaiHng models on proces- 
sive motor movement require the two motor domains to 
work in sequence, alternating attachment states allowing 
the motor to move forward without completely releasing 
from the filament [33, IH, 113, El ■ However, some motor 
proteins such as Myosin IXb have only a single motor 
domain and yet is still a processive motor |35l 139 . l4a | . 
In addition, numerous other experiments have been per- 
formed on single headed molecular motors, both natu- 
rally occurring and derived single-headed mutants. They 
show a continued ability to move and produce a force 
against an external load despite the obvious structural 
deficiencies with regards to the models based on two-head 
mechanisms ^ El El E El • 

Apart from these recent experimental findings, new 
models for motor proteins based on the idea of Brown- 



ian motors have been intensively investigated during the 
last decade HHH El Ullllj . When Brownian parti- 
cles are in a periodic but spatially asymmetric potential 
and away from thermal equilibrium, they can generate 
mechanical work by rectifying thermal fluctuation |52l| . 
Since motor proteins move along an asymmetric filament 
and ATP hydrolysis creates nonequilibrium conditions, 
it has been proposed that motor proteins could be de- 
scribed by a Brownian motor model [43, |4^ ]^ |50, Isil ■ 
A Brownian motor based on the flashing ratchet model 
moves as a step motor in consistent with experiments. 
However, Brownian motors also provide for the possi- 
bility of moving more than one step per chemical cycle 
in contrast to the powerstroke model. Furthermore, the 
step length of the Brownian motor is independent on the 
neck length. These features possibly account for some of 
the experimental data which the powerstroke model fails 
to explain |30ll53j . 

The powerstroke model and Brownian motor model 
has been thought of as two opposing extremes in the 
search for motility mechanisms of molecular motors. In 
the former, thermal fluctuations have no role at all. On 
the other hand, the latter is incapable of producing mo- 
tion without thermal fluctuations. As a result, it was 
commonly believed that only one of them is the correct 
mechanism. However, it has recently been postulated 
that both mechanisms may be involved in the actual mo- 
tor proteins [l^, |32| • 

The purpose of the present work is to develop a sin- 
gle mathematical model that includes both powerstroke 
and Brownian motor mechanisms. Since the coupling be- 
tween the rotational and linear translational degrees of 
freedom is essential to the powerstroke mechanism, our 
first step is to develop a model which couples rotational 
motion to linear motion (Section^. In the first model 
( Section ^nj, we demonstrate that free thermal diffusion 
in the rotational degree of freedom is sufficient to induce 
directed linear motion via a flashing ratchet mechanism. 
A model based on the behavior of motor proteins that is 
discussed in Section Hvl shows directed translational mo- 
tion driven by either a Brownian motor or a powerstroke 
mechanism depending on the parameters used. A more 
realistic motor protein model involving three chemical 
states is proposed in section |V] This motor is simultane- 
ously capable of both powerstroke and Brownian motor 
mechanisms during a single chemical cycle. Relevance 
of the present models to actual biological motors is dis- 
cussed in the final section. 



II. ROTATION-TRANSLATION COUPLING 

While various mathematical models for molecular mo- 
tors driven by rotation-translation coupling have been 
proposed in the past|5J,|53, these models are not closely 
related to known biological molecular motors. A more re- 
alistic model is necessary to explain the experimentally 
observed behaviors of biological molecular motors, such 




FIG. 1: Model for molecular motors with a coupled rotational 
and translational mechanism for producing directed motion. 
(a) Unbound state. Unattached of its filament track, the mo- 
tor is free to thermally diffuse in both x and 0. The resulting 
probability distribution A is shown in the lower image. In 
principle, Dx j^ Dg so A is depicted as an ellipse. The dashed 
lines depict the location of the potential barriers in the bound 
state for a point of reference, (b) Bound state. Upon binding 
to its filament, the pivot point P of the motor shifts towards 
the binding site, activating the coupling between rotation and 
translation. The potential contours shown in the lower image 
corresponding to the bound state, depict the probability of 
remaining in the same potential well at A or moving to the 
next at B. As can be seen here, an asymmetric potential U{x) 
is required to bias the motion of the motor. 



as step size distributions and velocity-force relationships. 
Here, we try to develop a mathematically simple, yet 
sufficiently realistic model by taking into account exper- 
imental observations as much as possible. 

Motor proteins undergo two major conformational 
changes during a chemical cycle, first when ATP hydrol- 
ysis takes place and second when ADP is released from 
the protein. The main change occurs between the head 
and neck of the protein such that the angle between the 
neck with respect to the head changes. For simplicity, we 
consider a simple rod of length ^ = 2r as a motor particle 
whose degrees of freedom are defined as the position of 
the center of mass x and a rotational angle around the 
center (measured from the normal to the filament). Since 
free rotation does not generate translational motion, it is 
therefore necessary to couple the rotational motion with 
the translational degree of freedom. Since the release of 
ADP occurs when the protein is attached to the filament, 
a large conformational change in the neck associated with 
the release would produce a rotation about the point of 
attachment to the filament. The resulting motion about 
that pivot point produces translational motion of the cen- 



ter of the motor that is proportional to r sin 9 as shown 
in Fig. n 

The potential energy of the motor in the bound state 
is determined by the location of the contact point xp and 
the orientation of the motor 9. In general, these two de- 
grees of freedom are weakly coupled but not in a way to 
produce a translational motion out of the swing motion. 
Therefore, we assume an uncoupled form of potential en- 
ergy: 



V{xp, 9) = Urot{e) + Utransixp) 



(1) 



where Urot represents a conformational change in the 
neck region of the protein and Utrans is a binding en- 
ergy between the motor and the filament. Since the fil- 
ament is periodic, the binding energy is also periodic as 
Utrans{xp + L) — U transit p) whcrc L is the pcriod of the 
filament. Noting that Xp and center-of-mass coordinate x 
are related by xp — x — r sin 9 k, x — r9 |5y| , the binding 
energy as a function of x and 9 is given by 



V{X, 9) = Urot{0) 4- Utrans{x - rb 



(2) 



This potential energy function provides the desirable 
form of coupling between rotational and translational de- 
grees of freedom, allowing for the conversion of rotational 
motion to translational motion as shown below. 



III. A TWO-STATE BROWNIAN RATCHET 
DRIVEN BY ROTATIONAL DIFFUSION 

Although motor proteins in general have several chem- 
ical states during ATP hydrolysis, we will restrict this 
discussion to a two-state model describing a bound and 
unbound state for a motor protein. When a potential is 
absent (unbound state) the motor freely diffuses in its 
two degrees of freedom, x and 9. As soon as the poten- 
tial is turned on (bound state), the motor drifts toward 
a local potential minimum. Switching between these two 
states will produce net movement in the x direction when 
the potential is asymmetric along x. To illustrate this 
ratchet mechanism, consider motors confined to an equi- 
librium point A of the potential well, as seen in Fig.^b). 
When the potential is removed, the motors diffuse over 
the shaded area in Fig. ^a). As the potential turns on 
again, a portion of the shaded area now belongs to the 
basin of next potential minimum B in Fig. H^b). The 
motors that reached this area advance to the new mini- 
mum B. The rest of the motors move back to the original 
minimum A. In the particular case shown in Fig. ^ the 
translational diffusion alone is too small to generate di- 
rected motion. However, the addition of the rotational 
diffusion brings the motors to the next basin. 

This two-dimensional Brownian ratchet mechanism 
can be put into a simple mathematical expression in a 
slow switching limit [48j. Assuming that the potential 
Utrans IS Sufficiently deep and the transitions between 
unbound and bound states are slow, the initial probabil- 
ity distribution of the motors is given by P(a;, 0, t = 0) = 



5{x — xa)5{9 — 9q). When the potential is removed, the 
distribution spreads with diffusion constants D^ and Dg 
as 



P{x,9,t) ^ Prot{e,t)Ptrans{x,t) 

1 / {x-xaf {9-9^Y 



2nty/DjJi 



: exp 



2D^t 2Dot 



(3) 



When the potential is reinstated at i = r, the probability 
that the motors leave the basin of the current potential 
wells in the forward (-fx) and backward (-x) direction are 
given by 



W+ir) 



oo p{x~d)/r 

dx / d9P{x,9,T) 

-oo J —OQ 



-erfc 



2 \^{2D^ + 2r^De)T ' 



(4) 



and 



nOQ n-\-OC 

W^{t) = I dx d9P{x,9,T) 

-oo J ~{x+L — d) /r 



= -erfc 



L-d 



2 \^{2D, + 2r^Dg)Tl ' 



(5) 



where d is the shortest distance between the potential 
minimum and maximum, [see Fig. 1(b)] and erfc is the 
complementary error function. If the motor is in the 
unbound state for period i„ and in the bound state for tt 
during one on-off cycle, then the average velocity of the 
motor is about 



tb + tu 



{W+{t^)^W-{tu)). 



(6) 



Here we neglected the possibility for the motors to diffuse 
beyond the nearest neighbor potential wells. When the 
potential is symmetric {d — L/2), W+ — W- and thus 
no net directed motion occurs. 

The above discussion can be further confirmed by a di- 
rect simulation of Langevin equations for the overdamped 
motor operated at a temperature T: 

1.^ = -^^^P^<^it) + V^k^W. Ut) + F (7a) 



dt 



dx 



7e 



d9 dV{x,9) r——- 

-jr^ T^ — (j{t) + V^kBTjg ^g{t), (7b) 



dt 



89 



where F is an external load and ^i is a Gaussian white 
noise define by 

UtKjit'))=S.,S{t-t') i,je{x,9). (8) 

The on-off function a{t) takes values and 1 alternatively 
with transition rates kon and koff- Although in general 
7a; would have an angular dependency according to the 
current orientation of the rod, here we assume that it 
takes constant values. 




FIG. 2: The relationship between average motor velocity 
and the motor radius r. The circle markers correspond to 
a— 0.0, the squares to a=1.0, and the triangles to q=10.0. 
The case a=0.0 corresponds to a purely translational ratchet. 
For nonzero a with a sufficiently large value for r, the ro- 
tational component of the motion provides a unique contri- 
bution to the motor velocity. The solid lines correspond to 
the long time limit velocity obtained from Eq.|S| Parameter 
values: t„=5.0, tt,=100.0, K=1.0, D=0.0075, d=0.36. 



As a simple model system we use a ratchet potential 
in the x direction: 



Utransix) = Uq 






1 . /47r 
— sm — X 
4 \ L 



(9) 



and a simple harmonic potential with a spring constant 
K for rotation: 



UrotiO) 



K 



(10) 



For simplicity, a value of K=l.{) will be used throughout 
the current study. 

To reduce the number of parameters, we normalize 
time, distance and energy as 



t = 



Upt 



X U_ 



(11) 



and accordingly parameters are normalized as f == ?'/L, 
K = K/Uq, F — FL/Uq. For simplicity, we will omit the 
tilde on the normalized quantities. Then, the equations 
of motion {Tj) become 



dx 

d0 
'di 



dV{x, d) 
dx 
dV{x,e) 



/2D ^^ + F 



a{t) + V2aD ^g 



(12a) 
(12b) 



where the diffusion constant is defined hy D = ksT/U^). 
A dimensionless constant a = r^^xHe is determined by 
the detailed shape of the motor. The contribution of ro- 
tational degree of freedom is controlled by the shape of 
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FIG. 3: The effect of applying a constant load F < to a mo- 
tor according to the equations of motion given by Eq. 1121 . 
(a) We see a distinct nonlinear behavior towards the applica- 
tion of an increasing external load F when the motion relies 
primarily upon the rotational ratchet mechanism. Parameter 
values: q:=10.0; -R'^1.0; kon=koff=0.5. (b) The nonlinear 
behavior disappears when equal proportions of both trans- 
lational and rotational ratchet are used. Parameter values: 
a=1.0; A'=1.0; kon=koff=0.5. 



the motor and hence the value of a. When a = 0, the 
rotational motion has no contribution and the present 
model is equivalent to a conventional flashing ratchet. 
For a > 1, the rotational diffusion plays a more signifi- 
cant role than the translational diffusion. 

Integrating the coupled Langevin equations in Eq. 112|l 
with the Heun method, we can explore the dynamics of 
this stochastic model. Examining the trajectories of nu- 
merous motors allows us to determine the average behav- 
ior of the system. For our purposes, an average over 500 
realizations of the system will be sufficiently accurate. 

Figure |5] shows that the velocity contains no depen- 
dence on r for a^O.O as expected. For nonzero values of 
a, the velocity begins to increase with the motor length r 
beyond the values expected from the translational ratchet 
alone, showing that the rotational component of the mo- 
tion introduces an additional mechanism for producing 
directed motion in the motor. The speed of the motor 
increases in a roughly linear fashion with respect to the 



motor length r similar to the way that most motor pro- 
teins are thought to behave. As the length of the motor 
increases, a change in the angle 9 will give the motor a 
longer reach, enabling it to move several potential peri- 
ods over one cycle with a sufficiently long motor. As seen 
in Fig. 12 the rotational motion contributes to the mo- 
tion with increasing r until the motor is long enough to 
reach the nearest neighboring potential well. Beyond a 
certain length, the velocity is expected to saturate, until 
it is then long enough to reach the next nearest neighbor. 
This behavior is corroborated by the expected behavior 
of the long time limit case for these various values of a. 
Figure 121 shows close agreement between Eq. |H1 described 
above and the stochastic simulations. 

By applying a constant load _F < 0, we can investigate 
the ability of this model to do work against a load. Figure 
O shows the results of two cases: a^lO.O where the ro- 
tational diffusion dominates and a=1.0 where rotational 
and translational diffusion equally contribute to the mo- 
tion. The limit where the motor moves primarily under 
rotational diffusion exhibits a nonlinear response to an 
applied load. Upon moving towards the limit of equal 
contributions {i.e. a—l.Q) the nonlinearity disappears. 
The nonlinearity is therefore due to the anisotropic dif- 
fusion. 



IV. A TWO-STATE BROWNIAN MOTOR WITH 
POWERSTROKE 

In the previous model, coupling between the rotational 
and translational degree of freedom takes place only when 
the motor protein is bound to an active site on the fil- 
ament. The rotation toward an equilibrium angle ad- 
vances the motor to a minimum within the basin of the 
potential well. We shall call this rotation the forward 
stroke. The forward stroke however cannot drive the 
motor over the potential barrier and only diffusion, ro- 
tational or translational, brings the motor to the next 
potential well. On the contrary, in the standard power- 
stroke model conformational change of the motor protein 
drives the motor without diffusion. The key point of the 
powerstroke model is the presence of the recovery stroke 
which the motor performs when it is unbound. Molecu- 
lar motors typically consist of a motor domain attached 
to a neck and tail region. The tail region typically binds 
to some sort of cellular cargo or, in the case of skele- 
tal muscle, the tail binds together with other proteins 
to form a bundle that works in tandem. As opposed to 
our free floating motor presented in the previous section, 
the motor protein does not freely rotate in the unbound 
state and instead is driven to an equilibrium angle by the 
recovery stroke. 

As an alternate model to our previous Brownian 
ratchet, we consider a motor which is not entirely free 
to move, but is tethered by one end to some external 
object assumed to be much larger than itself. When the 
motor is unbound, it rotates about the pivot point P' 




FIG. 4: A model for motor proteins with one end of a motor 
of length l tethered to an external cargo, (a) Unbound state. 
As before, the motor is free to diffuse in x and 6, but 6 drifts 
toward a stable angle 6\. Overlaid on the potential contours 
is the new probability distribution at B centered about 6i, 
away from the potential minimum of the bound state (shown 
at A'), (b) Bound state. On binding, the pivot P shifts, giving 
rise to the coupling term in Utransix). Here the stable angle 
for the bound state 62 is not necessarily equal to 9\. Given 
a large enough change in the protein conformation between 
the bound and unbound states t.t^Q=li(Q-2 — di), motion can 
occur deterministically without the need for the probability 
distribution at B. 



toward an equilibrium angle 61 (recovery stroke) as il- 
lustrated in in Fig. EJa). When the motor binds to the 
filament, the pivot point shifts to the other end of the 
motor, P in Fig.^b). The motor rotates to a new equi- 
librium angle 62 about this point (forward stroke). The 
two different equilibrium angles 9i and 62 represent con- 
formational changes in the motor proteins due to ATP 
hydrolysis. 

Suppose that the motor is located at a bound state A 
in Fig. 01b). Upon detaching from the filament, the re- 
covery stroke takes the motor from A' to B in Fig. 0Ja). 
If the conformational change, iA9 — £{62 — Ox) is large 
enough, the point B is in the basin of the next potential 
well. When it binds to the filament, the forward stroke 
bring the motor to C, generating a net translational mo- 
tion from A to C . This process is entirely deterministic 
and the directed motion is generated simply by alter- 
nating forward and recovery rotations without diffusion. 
However, the motor is still subject to large thermal fluc- 
tuation, which may hamper the powerstroke. Indeed, the 
thermal fluctuation brings some motors back to the basin 
of the original potential well as illustrated in Fig. ^ and 
thus reduces the average speed. On the other hand, the 
thermal fluctuations also make a positive contribution 
to the directed motion especially when the powerstroke 
fails. For example, one could imagine a motor prema- 



turely binding to the filament before the recovery stroke 
completes. In a deterministic model that powerstroke 
would fail to advance the motor to the next potential 
well. With the addition of the thermal diffusion, some 
of the motors that would fail the powerstroke are able to 
reach the next potential well with help of that diffusion. 
Furthermore, when the conformational change £A0 is too 
small to reach the next minimum, the Brownian motion 
is still capable of producing directed motion on average. 
Unlike the previous model, rotational diffusion in the 
unbound state is accompanied by a drift due to the con- 
formational change tA6. This rotational drift determines 
the direction of motion. When the sign of £A0 is reversed, 
the motor moves in the opposite direction. Recent exper- 
iments have exhibited such flux reversal for motors mod- 
ified to swing their neck in the opposite direction to the 
original motors |3J,|33. Since the symmetry is broken by 
these conformational changes, it is not necessary for the 
potential to break the symmetry along the x axis. A sim- 
ple periodic potential in the normalized units described 
earlier will suffice for our purposes: 



Utrans{x) = COs(27rx). 



(13) 



The net potential for the unbound and bound states are 
respectively given by 



Vi{x,e) 



^( 

2 ^ 

^2(X,0) = |( 



hf + Utrans{x-ie), 



(14a) 
(14b) 



where x is now the position defined by the neck of the 
motor protein at P' . The movement of the motors is de- 
termined by the Langevin equations in normalized units: 



dx_ _ dV{x,e,t) 

dt dx 

d9 dV(x,0,t) fTT-R . 



(15a) 
(15b) 



where the time-dependent potential is defined as 
V{x, 0,t) = yi(x, e)[l - a{t)] + V2{x, e)a{t). 

We begin by examining the relative change in the con- 
formation iA9 between the two states Vi and V2. Fig- 
ure [SI shows the average velocity of the motor as £A9 
is increased for several different values of the diffusion 
constant D. In the case of D=0.0, the thermal fluctua- 
tions are disabled. For a small neck length f or a small 
A9, no velocity is observed. However, once the critical 
value of £A9 w 0.5 is reached, the motor begins to move 
without the need of the thermal fluctuations. This de- 
terministic motion corresponds to a powerstokc motion, 
related both to the neck length of the motor given by £ 
and the angle through which the motor rotates between 
the unbound and bound states. Although in principle 
the motor moves forward one potential through each hy- 
drolysis cycle, close to the critical value of £A6 the motor 
would require a very long switching time to relax to the 
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FIG. 5: (a) The average velocity of the motor (x) as a func- 
tion of the change in the stable conformations £A9 between 
the two states with the unbound angle 9i fixed. Parameter 
values: tu=tt^4.0; K=1.0; 6'i=0.0; q=1.0. (b) The average 
velocity as a function of A6 with the bound stable angle 82 
fixed. Parameter values: t„=ii,=4.0; 7^=1.0; D=0.01; a=1.0. 



potential minimum of the bound state in order to achieve 
that expected velocity. 

In the case of £> ^ 0, Fig. [S] shows that even for small 
values of D motion can be seen in the region of £A9 < 0.5 
where the deterministic powerstroke was not possible. 
The thermal fluctuations of the medium surrounding the 
motor can be rectified to produce a nonzero average ve- 
locity. The mechanism for rectifying the thermal fiuctu- 
ation follows the Brownian ratchet mechanism from the 
coupled rotational and translational coupling described 
in Section imi 

Also of interest is the case where £A9 < 0. This pa- 
rameter regime would be equivalent to a motor taking 
a power stroke in a negative direction, and as a results 
Fig. El shows an average negative velocity for the motor. 
This observation is of potential interest in describing the 
motion of motor proteins such as Myosin VI, observed to 
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FIG. 6: (a) An analysis of the tether motor system showing 
both the case of smaU and large values of D. The markers cor- 
respond to the Langevin simulation and the lines correspond 
to the long time limit theory. Parameter values: fu=ib=16.0; 
iS'=1.0; a=1.0. (b) An external load J^ < applied to a 2 
state tethered motor system. Parameter values: fu=i(,=16.0; 
J^=1.0; D=Q.01; a=1.0. 



t/S.9 becomes very large (£A0 > 1.0), the determinis- 
tic motion dominates the model. For small changes in 
the angle {t/S.9 < 0.5), the Brownian ratchet mechanism 
dominates. Between these two limits a combination of 
the two mechanisms contribute to the average velocity 
of the motor. The resulting model forms a very robust 
motor that can produce directed motion under a wide 
variety of parameters. 

As before, the ratcheting motion found in this model 
can be studied analytically in the long time limit as de- 
scribed before with little difference. However the limits 
of the integral in Eq. Q change to take into account 
the change the angle 9i between states. The new forms 
for W+ and W- can be obtained by replacing d with 
L/2 - iA9 in Eqs. ^ and ^. A nonzero iA9 breaks 
the symmetry of the system and gives rise to W-f ^ W- 
even though we are using a symmetric potential. The 
resulting current then follows the same form as stated in 
Eq. ((BJ. In comparing the results of the long time limit 
calculation to the stochastic simulations in Fig. IHIa), we 
again find close agreement to the motion seen in the long 
time limit. 

Motor proteins typically bind a large cargo to their 
tail in order to transport it from one part of a cell to 
another. As before, we can apply an external load F < 
to these motors to study their ability to pull such a load 
that might be attached to it. In Fig. Elb), the same 
type of linear behavior is observed as before with the 
previous model. However, Figure El^b) also shows that 
as the value of iA9 is increased the maximum force with 
which the motor pulls increases. It has been postulated 
that a Brownian ratchet is unlikely to be a mechanism for 
motor proteins as those models tend to produce forces too 
small to account for the forces observed experimentally 
[l|. Based on the present results, it appears that as we 
increase the contribution of powerstroke to the function 
of the motor, the ability of the motor to transport cargo 
more effectively also increases. 



A THREE STATE MODEL 



move in the opposite direction as the plus-ended Myosin 
V proteins. Using the angle 6*1=0.0 for the unbound con- 
figuration gives a completely symmetric response given 
the symmetric nature of Utrans{x). On the other hand, 
by fixing the bound stable angle 02 to some non-zero 
value, the symmetry disappear. To meet the conditions 
for iA9 < with 92 fixed, 9i must become increasingly 
large to obtain any negative motion. The harmonic po- 
tential term in Eq. H14I) quickly becomes too large to 
sustain negative motion and drives the motor back to- 
wards the positive. It is likely that this case is not well 
represented by the current model and this behavior is an 
artifact of the approximations that Urot is a harmonic 
potenial and rsa\9 w r9. 

In summary, the diffusionless curve in Fig. |S1 shows 
the contribution of the deterministic motion. When the 



The cyclic chemical reaction that myosin undergoes 
as it hydrolyzes ATP can be expressed with varying de- 
grees ofprecision involving numerous substeps of the re- 
action 0, ISTl l58| . Detailed analysis of the structure of 
myosin SI has indicated the existence of at least three 
distinct conformational states dependent on the state of 
the bound ATP molecule. H S Q, H i, El . In addition, 
some myosin working strokes have been resolved into at 
least two distinct substeps [2^ jSj, |63 , yielding more ev- 
idence of multiple conformations of the motor protein. 
Here, we will utilize three distinct conformational states 
with a chemical reaction scheme similar to the simple 
mechano chemical trigger model put forth by Keller and 
Bustamante [53 . This setup enables us to create a hy- 
brid motor that contains elements of both a Brownian 
motor as well as a deterministic powerstroke. 
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FIG. 7; A three state reaction sdieme for a motor protein 
and a typical trajectory for such a motor. Rate constants gov- 
erning the reaction are labeled as kij and the conformational 
angle for each state is 9i where i,j € (1,2,3). Trajectory pa- 
rameter values: D=0.01; a=1.0; fci2=0.5, fc2i=0.2, ^23=0. 2, 
fc3i=0.004; (72=0.2, {/3 = 1.0; K^l.O; ^6*12=0.2, M6l23 = 1.0. 



Following the cyclic reaction scheme shown in Fig. \7\ 
the transition between a motor containing ATP in an 
unbound state and its subsequent weakly bound state to 
an actin filament as the ATP is broken down into ADP 
and Pi is reversible. This transition induces a small con- 
formational change iA9i2 = £{62 — Oi) which is not big 
enough to perform a powerstroke. However a back and 
forth transition between this unbound and weakly bound 
state can generate directed motion by rectifying rota- 
tional diffusion as discussed in the previous sections. The 
transition from a weakly bound to a strongly bound state 
involves the irreversible release of the inorganic phos- 
phate Pi and is accompanied by a larger change in the 
conformational state £A923 — ^(^3 — ^2) providing the 
powerstroke. Completing the cycle involves the release of 
ADP and binding of a new ATP molecule, returning the 
motor protein to an unbound state. Cycling along these 
three states, the motor has the capability of utilizing both 
the Brownian motor and powerstroke mechanisms. 

Comparing the typical trajectory of this model shown 
in Fig. 13 to single molecule experiments performed in 
recent years shows a possible mechanism explaining ob- 
served behavior. The model presented here exhibits oc- 
casional backward steps and multiple forward steps dur- 
ing the Brownian ratchet phase similar to the events ob- 
served by Yanagida's group |23,|30|. The additional pow- 
erstroke as the motor enters the strongly bound state 
provides another step forward along the potential, sim- 
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FIG. 8: (a) Motor velocity as a function of varying confor- 
mational changes from state 1 to 2, and from state 2 to 3. Pa- 
rameter values: D=0.01; a— 1.0; fci2=0.5, ^21=0.2, ^23=0. 2, 
itsi =0.004; C/2=0.2, t/3=1.0; 7^=1.0. (b) Motor velocity as a 
function of an applied load -F < 0. Parameter values: as in 
(a) with M6l23 = 1.0. 



ilar to the two stag e working stroke observed in some 
experiments |26ll32.l60| . 

Figure |HIa) demonstrates the change in velocity as 
both the magnitude of the powerstroke and the confor- 
mational change associated with the Brownian motor are 
increased. In the case of iA6i2=0, the Brownian motor 
mechanism makes no contribution to the velocity of the 
motor. As iA6i2 is increased, the Brownian motor mech- 
anism has the capability of tripling the overall speed of 
a motor, implying that the motor could take on aver- 
age two extra steps during the Brownian motor phase of 
the hydrolysis cycle with adequate tuning of the parame- 
ters. As seen in Fig.lHJb), this three state system has the 
added advantage of withstanding a larger external load 
opposing the motion than the simpler two state model. 



VI. DISCUSSION 

Conventional Brownian ratchets encounter various dif- 
ficulties when applied to biological motors. In the present 
model, coupling the translational and rotational de- 
grees of freedom preserves the conventional translational 
ratchet mechanism, but also provides a means to produce 
translational motion through rectifying thermal fluctua- 
tions in the rotational degree of freedom. Further, with 
the addition of an asymmetric conformational change be- 
tween bound and unbound states, we obtain a model ca- 
pable of rectifying thermal fluctuations without the need 
for an asymmetric potential. In the conventional Brow- 
nian ratchet model the direction of the motion is deter- 
mined by a spatial asymmetry in the translational degree 
of freedom. However, in the present model, the direction 
of motion is determined by the asymmetry in the confor- 
mational change in the motor proteins rather than the 
asymmetry in the filament, consistent with experimen- 
tal observation |34L l37l |. If this conformational change 
between the bound and unbound states is large enough, 
the motor gains the added ability to move by powerstroke 
without diffusion. Depending on the parameters chosen, 
either mechanism can be utilized ultimately gaining the 
benefits of both. 

With the minimal three chemical state description de- 
tailed in the previous section, we introduce a realistic 
model for single headed motor proteins, which captures 
various qualitative properties of motor proteins. Recent 
experiments on monomeric Myosin I have emphasized the 
importance of both the neck length and the degree of 
neck rotation in determining the step size of a motor pro- 
tein [13. The present model takes both factors into ac- 
count in defining the coordinate of conformational change 
as IM). Other single molecule experiments have resolved 
at least two distinct conformational changes during the 
working stroke of Myosin I 60] , one corresponding to a 
displacement of '-^ 6.5 nm and another of ~ 5.5 nm, both 
very close to the size of the actin filament repeat. The 
multiple steps observed in Myosin I can be explained by 
the different phases of motion through an ATPase cy- 
cle of the present model. As seen in the trajectory in 
Fig. [7| a diffusive step followed by a powerstroke would 
each produce displacement in a series of steps equal to 
the filament periodicity L. 

Many single molecule experiments with Myosin II have 
observed step sizes on the order of 4 — 10 nm [20, yfl 
l6ll l6a, |63, ISJ, |63. However, according to Tanaka et 
al. "63] , a random orientation of myosin with respect to 
an actin filament produces the observed average working 
stroke of '-^ b.5nm. When myosin is close to its natural 
orientation with respect to actin as it would be in a whole 
fiber, the step size is closer to three times that value, 
or ~ 15 nm. Kitamura et al. |30j observed substeps of 
~ 5.3 nm along actin subunits. As the dwell time of 
these substeps had no dependence on ATP concentration, 
it was proposed that an average of 2.5 substeps during 
a single ATPase cycle makes up the observed '^ 13 nm 



step per ATP molecule consumed. A simple deterministic 
motor protein should take a single step per ATP molecule 
consumed. However, as illustrated in Fig. ISfa), a motor 
protein with £A6'i2=0.5 on average should take roughly 
two 5-nm diffusive steps along the filament subunits in 
addition to the powerstroke during the consumption of a 
single ATP for a total step size of 15 nm. 

In the case of dimeric, processive motor proteins like 
Myosin V, models for the processivity involve some type 
of coordination between the two motor domains of the 
myosin. Investigations into single-headed Myosin V mo- 
tors have revealed a directional force de pen dence on the 
kinetics that regulate ATP hydrolysis psl l4a |. These 
force dependent kinetics indicate a possible method of 
communication between the two individual motor do- 
mains. When the leading head of the motor pair finds its 
binding site on the actin filament and enters a strongly 
bound state, the strain caused on the trailing head en- 
courages its release from the filament. At the same time 
the release of the leading head is slowed, effectively an- 
choring the motor protein to the filament allowing the 
trailing head to swing forward and find its next binding 
site. While this motion of the trailing head forward is 
thought as a powerstroke process [ifl IM, H^ ll^j El, 
some recent experiments suggest an alternative. The 
working stroke of a single-headed Myosin V motor pro- 
tein produces a displacement of only ~ 25 nm, signifi- 
cantly shorter than the observed step size of 36 nm for 
wild type two-headed Myosin V |33| . The motion of this 
wild type Myosin V has since been resolved into a 12 nm 
step followed by a 24 nm step [63 • Due to the insuf- 
ficient length of the working stroke of a single-headed 
Myosin V motor, it has been proposed that the power- 
stroke of the leading head partially positions the trailing 
head and the remaining distance is covered by a diffusive 
process 11, 32] . Others have shown that single headed 
Myosin V molecules can move processively with the mo- 
tor protein held in close proximity to the filament |4j|. 
In this study it was proposed that through a series of dif- 
fusive steps utilizing a strain sensors similar to the one 
demonstrated by Veigel et al. j4^ and Purcell et al. [4g, 
a single motor protein can achieve a step size of 32 nm. 

Myosin VI, another dimeric motor protein, has a neck 
length much less than that of Myosin V, yet still produces 
a step length of similar size j2g,^l|. As a result, any de- 
terministic powerstroke of this motor protein probably 
has little to do with its step and is most likely driven 
by thermal motion. The present model provides a single 
mechanism allowing a motor to take advantage of both 
diffusion as well as a subsequent powerstroke/anchoring. 
Given the large neck of Myosin V, a significant power- 
stroke on the order of 25 nm would be reasonable based 
on a large value for i, followed by a Brownian motor 
phase positioning the trailing head at its full step distance 
from the leading head. On the other hand, a model for 
Myosin VI could rely heavily upon the Brownian motor 
mechanism in order to reproduce the observed behav- 
ior. We will further explore these dimeric motor proteins 
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based on our model in a future paper |67j . 

Some monomeric motor proteins such as Myosin IX 
are thought to involve a diffusive step in order to achieve 
some measure of processivity jSa, |4Z, 68, o!^. However, 
to maintain this processivity monomeric motor proteins 
must remain in the vicinity of a filament while in a dif- 
fusive state. Both Nalavadi et al. ^8] and Kambara et 
al. |69| have proposed a mechanism. The structure of 
Myosin IX contains a unique insertion in the motor do- 
main that may harbor a second binding site. This in- 
sertion could tether the protein to its filament while still 
allowing enough freedom for the protein to diffuse. The 
motor protein then diffuses towards its next binding site 
before entering another powerstroke. Given the assump- 
tion of this special tether in Myosin IX, the present model 
provides a mechanism to efficiently guide the motor to its 
next binding site by rectifying thermal fluctuations be- 
fore entering a strongly bound state. 

All of these single molecule constructs have been 
observed to pro duce a force on the order of a few 
pN [llillillie'.'ei'.'ellTOll. Examining the normalized 
load F from Ea. (|12|l . for L=5A nm, the approximate 
length of an actin subunit, and Uo=20kBT one unit of 
normalized force is about 15.8 pA^. The maximum forces 



generated by the monomers shown in Fig. |HIb) are then 
approximately 2.6 — 3.5 piV, comparable to the forces ob- 
served experimentally. 

In conclusion, we have developed a realistic model of 
motor proteins which simultaneously includes the tra- 
ditional Brownian ratchet model and the powerstroke 
model in a simple unified form. The properties of this 
model are qualitatively consistent with recent experimen- 
tal results. Quantitative investigation is currently under- 
way. Furthermore, we are applying the present model to 
dimeric processive motors as well as muscle fiber. Al- 
though the model is developed for the myosin superfam- 
ily, we believe it can be applied to other molecular motors 
such as kinesin. 



VII. ACKNOWLEDGMENTS 

The authors would like to acknowledge Kimmie Farris 
and Erin Darnell for their help in this work. This work 
was partially supported by National Science Foundation 
(NSF) GK-12 Award No. 0139108. 



[1] J. Howard, Mechanics of Motor Proteins and the Cy- 

toskeleton (Sinauer Associates, Inc., Sunderland, MA, [18 
2001). 
[2] A. F. Huxley, Prog. Biophys. Biophys. Chem. 7, 255 [19 

(1957). 
[3] A. F. Huxley and R. M. Simmons, Nature 233, 533 [20 

(1971). 
[4] M. J. Tyska and D. M. Warshaw, Cell Motility and the [21 

Cytoskeleton 51, 1 (2002). 
[5] I. Rayment, W. R. Rypniewski, K. Schmidt-Base, et al., [22 

Science 261, 50 (1993). 
[6] M. Whittaker, E. M. Wilson-Kubalek, J. E. Smith, et al., [23 

Nature 378, 748 (1995). 
[7] J. D. Jontes, E. M. Wilson-Kubalek, and R. A. Milligan, [24 

Nature 378, 751 (1995). 
[8] A. Houdusse, V. N. Kalabokis, D. Himmel, A. G. Szent- [25 

Gyorgyi, and C. Cohen, Cell 97, 459 (1999). 
[9] A. Houdusse, A. G. Szent- Gyorgyi, and C. Cohen, PNAS [26 
97, 11238 (2000). 
[10] J. N. Forkey, M. E. Quinlan, M. A. Shaw, J. E. T. Corrie, [27 

and Y. E. Goldman, Nature 422, 399 (2003). 
[11] M. A. Geeves, Nature 415, 129 (2002). 
[12] M. Anson, M. A. Geeves, S. E. Kurzawa, et al., EMBO 

J. 15, 6069 (1996). [28 

[13] T. Q. P. Uyeda, P. D. Abramson, and J. A. Spudich, 

PNAS 93, 4459 (1996). [29 

[14] T. Sakamoto, F. Wang, S. Schmitz, et al., J. Bio. Chem. 

278, 29201 (2003). 
[15] T. J. Purcell, C. Morris, J. A. Spudich, et al., PNAS 99, [30 

14159 (2002). 
[16] D. H. Schott, R. N. CoUins, and A. Bretscher, J. Cell. [31 

Biol. 156, 35 (2002). 
[17] C. Ruff, M. Furch, B. Brenner, et al., Nat. Struc. Biol. [32 



8, 226 (2001). 

D. M. Warshaw, W. H. Guilford, Y. Freyzon, et al., J. 

Biol. Chem. 275, 37167 (2000). 

D. Kohler, C. Ruff, E. Meyhofer, and M. Bahler, J. Cell 
Biol. 161, 237 (2003). 

A. Ishijima, H. Kojima, T. Funatsu, et al.. Cell 92, 161 

(1998). 

S. Itakura, H. Yamakawa, Y. Y. Toyoshima, et al., 

Biochem. Biophys. Res. Commun. 196, 1504 (1993). 

K. M. Trybus, E. Krementsova, and Y. Freyzon, J. Biol. 

Chem. 274, 27448 (1999). 

C. Perreault-Micale, A. D. Shushan, and L. M. Coluccio, 

J. Biol. Chem. 275, 21618 (2000). 

K. Homma, J. Saito, R. Ikebe, et al., J. Biol. Chem. 275, 

34766 (2000). 

H. Tanaka, K. Homma, A. H. Iwane, et al.. Nature 415, 

192 (2002). 

I. Lister, S. Schmitz, M. Walker, et al., EMBO J. 23, 

1729 (2004). 

S. Nishikawa, K. Homma, Y. Komori, M. Iwaki, 

T. Wazawa, A. Hikikoshi Iwone, J. Saito, R. Ikebe, 

E. Katayama, T. Yanagida, et al., Biochem. Biophys. 
Res. Commun. 290, 311 (2002). 

R. S. Rock, S. E. Rice, A. L. Wells, T. J. Purcell, J. A. 

Spudich, and H. L. Sweeney, PNAS 98, 13655 (2001). 

A. H. Iwane, H. Tanaka, S. Morimoto, A. Ishijima, and 

T. Yanagida, Journal of Molecular Biology 353, 213 

(2005). 

K. Kitamura, M. Tokunaga, A. H. Iwane, et al.. Nature 

397, 129 (1999). 

M. Y. All, K. Homma, A. H. Iwane, et al., Biophys. J. 

86, 3804 (2004). 

C. Veigel, F. Wang, M. L. Bartoo, J. R. Sellers, and J. E. 



11 



MoUoy, Nat Cell Biol 4, 59 (2002), ISSN 1465-7392. [52 

[33] R. D. Vale, J. Cell Biol. 163, 445 (2003). [53 

[34] A. L. Wells, A. W. Lin, L.-Q. Chen, D. Safer, S. M. Cain, 

T. Hasson, B. O. Carragher, R. A. Milligan, and H. L. [54 

Sweeney, Nature 401, 505 (1999). [55 

[35] A. Inoue, J. Saito, R. Ikcbe, and M. Ikebe, Nat. Cell Biol. [56 

4, 302 (2002). 

[36] C. B. O'Connell and M. S. Mooseker, Nat. Cell Biol. 5, 

171 (2003). [57 

[37] G. Tsiavaliaris, S. Fujita-Becker, and D. J. Manstein, Na- [58 

ture 427, 558 (2004). [59 

[38] S. Leibler and D. A. Huse, J. Cell Bio. 121, 1357 (1993). 
[39] A. Yildiz, J. N. Forkey, S. A. McKinney, et al.. Science [60 

300, 2061 (2003). 
[40] M. J. Tyska and M. S. Mooseker, Trends Cell Biol. 13, [61 

447 (2003). 
[41] A. Vilfan, Biophys. J. 88, 3792 (2005). [62 

[42] P. L. Post, M. J. Tyska, C. B. O'Connel, K. Johung, 

A. Hayward, and M. S. Mooseker, J. Biol. Chem. 277, [63 

11679 (2002). 
[43] T. Yanagida, K. Kitamura, H. Tanaka, A. H. Iwane, and [64 

5. Esaki, Curr. Opin. Cell Biol. 12, 20 (2002). 

[44] T. M. Watanabe, H. Tanaka, A. H. Iwane, S. Maki- 

Yonekura, K. Homma, A. Inoue, R. Ikebe, T. Yanagida, [65 

and M. Ikebe, PNAS 101, 9630 (2004). 
[45] C. Veigel, S. Schmitz, F. Wang, and J. R. Sellers, Nat 

Cell Biol 7, 861 (2005). [66 

[46] T. J. Purcell, H. L. Sweeney, and J. A. Spudich, PNAS 

102, 13873 (2005). [67 

[47] M. O. Magnasco, Phys. Rev. Lett. 71, 1477 (1993). [68 

[48] R. D. Astumian and M. Bier, Phys Rev Lett 72, 1766 

(1994). [69 

[49] R. D. Astumian, Science 276, 917 (1997). 
[50] F. Jiilicher, A. Ajdari, and J. Prost, Rev. Mod. Phys. 69, [70 

1269 (1997). 
[51] F. Jiilicher and J. Prost, Phys Rev Lett 75, 2618 (1995). 



P. Reimann, Phys. Rep. 361, 57 (2002). 

T. P. Terada, M. Sasai, and T. Yomo, PNAS 99, 9202 

(2002). 

Y.-X. Li, Physica A 251, 382 (1998). 

M. Porto, Phys. Rev. E. 63, 030102(R) (2001). 

Here we assumed 9 « 1 just for symplicity. Although 

actual angle may not be small, the main conclusions in 

the present work are valid without this approximation. 

D. Keller and C. Bustamante, Biophys. J. 78, 541 (2000). 

R. D. Astumian, Appl. Phys. A 75, 193 (2002). 

P. D. Coureux, H. L. Sweeney, and A. Houdusse, EMBO 

J. 23, 4527 (2004). 

C. Veigel, L. M. Coluccio, J. D. Jontes, et al., Nature 

398, 530 (1999). 

J. E. MoUoy, J. E. Burns, J. Kendrick- Jones, R. T. 

Tregear, and D. C. S. White, Nature 378, 209 (1995). 

A. D. Mehta, J. T. Finer, and J. A. Spudich, PNAS 94, 
7927 (1997). 

H. Tanaka, A. Ishijima, M. Honda, K. Saito, and 

T. Yanagida, Biophys. J. 75, 1886 (1998). 

M. J. Tyska, D. E. Dupuis, W. H. Guilford, J. B. Pat- 

lak, G. S. Waller, K. M. Trybus, D. M. Warshaw, and 

S. Lowey, PNAS 96, 4402 (1999). 

N. M. Kad, A. S. Rovner, P. M. Fagnant, P. B. Joel, 

G. G. Kennedy, J. B. Patlak, D. M. Warshaw, and K. M. 

Trybus, J. CeU Biol. 162, 481 (2003). 

S. Uemura, H. Higuchi, A. O. Olivares, E. M. De La Cruz, 

and S. Ishiwata, Nat. Struc. Mol. Biol. 11, 877 (2004). 

B. Geislinger and R. Kawai, in preparation. 

V. Nalavadi, M. Nyitrai, C. Bertolini, N. Adamek, M. A. 

Gocves, and M. Bahler, J. Biol. Chem. 280, 38957 (2005). 

T. Kambara and M. Ikebe, J. Biol. Chem. 281, 4949 

(2006). 

A. D. Mehta, R. S. Rock, M. Rief, J. A. Spudich, M. S. 

Mooseker, and R. E. Cheney, Nature 400, 590 (1999). 



